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Abstract 

In two previous papers, the Kerman-Klein-Donau-Frauendorf (KKDF) model 
was used to study rotational bands of odd deformed nuclei. Here we describe 
backbending for odd nuclei using the same model. The backbending in the 
neighboring even nuclei is described by a phenomenological two band model, 
and this core is then coupled to a large single-particle space, as in our previous 
work. The results obtained for energies and Ml transition rates are compared 
with experimental data for ^^^Lu and for energies alone to the experimental 
data for ^''^W. For the case of ^^^Lu comparison is also made with previous 
theoretical work. 
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I. INTRODUCTION 



In two previous applications the Kerman-Klein- Donau-Frauendorf (KKDF) model 
was used to study rotational bands of selected odd deformed nuclei. In both applications 
the system was described by an effective interaction which includes a monopole pairing and 
a quadrupole-quadrupole interaction. In the first application a large single-particle space 
was coupled to the ground-state band of the neighboring nuclei, whereas in the second 
application the same large single-particle space was not only coupled to the ground-state 
band but also to some of the excited bands of the core. 

As a second class of applications, we want to describe backbending for odd nuclei using 
the KKDF model. As has been known for more than two decades, deformed nuclei commonly 
show a rotational anomaly (known as backbending), where the energetically favored, or 
yrast, collective band undergoes an abrupt increase in its moment of inertia (as a function 
of frequency, for example). The generally accepted interpretation is that backbending of 
an even nucleus occurs when two neutrons (or protons) with high j break their pairing 
bond and rotationally align perpendicular to the symmetry axis 0. The rotational-aligned 
sequence of states is called the s-band. In this application, instead of following a microscopic 
approach to backbending that incorporates this physics, we have chosen to take a purely 
phenomenological approach using a model of two coupled bands. These will then be coupled 
to the extra odd particle in order to achieve a description of the corresponding phenomena 
in the neighboring odd nuclei. The reason for utilizing a phenomenological description of 
the even cores is that this provides a more accurate fit to the data than existing microscopic 
calculations. 

In the next section we will present the description of the cores, including the phenomeno- 
logical model and the results. In Sec. |IT| the results of the energy calculations for ^^^W will 
be presented, and in Sec. ^ the results for both energies and Ml transitions will be given 
for ^^^Lu. There is a brief concluding section. 

II. CORE PHENOMENOLOGY 

In this section we develop a phenomenological description of the backbending phe- 
nomenon for the even neighboring nuclei. It is not our purpose to build a complete and 
sophisticated model but to reproduce the energy levels and the values of the quadrupole 
matrix elements as accurately as possible, because the results of these calculations will be 
used as input to the calculations for the odd nuclei. The purpose of the model is simply to 
provide input data to our theory when experimental values are not sufficiently abundant. 

Consider two bands, the ground-state band and an excited band that "cross" at a certain 
angular momentum /. From experimental observations we know that if such a situation 
occurs, the bands repel each other (avoided crossing). From the mixture of the bands, we 
can conclude that the Hamiltonian of the system is in general not axially symmetric, and 
as a result, the projection K is not a good quantum number. We write the Hamiltonian of 
such a system as 

H = Ho + Hi, (1) 
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where Hq is the axially symmetric part of Hamiltonian and Hi is the K quantum number 
non-conserving Hamiltonian which gives rise to the avoided crossing. The unperturbed 
Hamiltonian, Hq, describes the ground-state rotational band and one excited band. The 
excitations in such bands can have the form of the simple rigid rotor, or of the variable 
moment of inertia model (VMI model) 0,^, or even more general form (see below). For the 
uncoupled systems, we take the angular momentum /, its projection M, and the projection 
on the body fixed axis K as good quantum numbers. Therefore, we denote the eigenstates of 
the unperturbated system as \IMK). Expressed in the unperturbated basis, the Hamiltonian 
for a given angular momentum will take the following form. 
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where I = I{I + 1) and Eg is the bandhead energy. The perturbation has, for the mo- 
ment, a general / dependence. The functions /«(/) which describe the angular momentum 
dependence of the uncoupled band i {i = 1, 2) are often represented as a polynomial in 
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(3) 



with expansion coefficients X°, X^. This conventional expansion has poor convergence prop- 
erties at large I. A better expansion of energy as a function of angular momentum is given 
by the variable moment of inertia model (VMI) in which the moment of inertia is a 
function of angular momentum 



(4) 



with two parameters C and Xq. The variable moment of inertia X(/) is determined through 
use of the variational condition 



dE{I) 



dl{l) 



0. 



(5) 



This model can be shown to be equivalent to the two-term approximation to the Harris 
(cranking) formula ||^. The latter is an expansion of the energy in powers of the "rotational 
frequency" u, given by 



E{u;)=A'>u;^ + A'u;^ + A'u;'^ + 



(6) 



where , A^ etc are parameters that are fitted to experimental values. Since the concept 
of continuous angular frequency is not well defined for quantum mechanical systems, the 
definition used is derived from the corresponding classical definition. In classical mechanics 
we have. 



dE{I) 
dl ' 



(7) 
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where E{I) is the energy at angular momentum /. In quantum mechanics we take a cue 
from this definition and choose the discrete definition 



where J+i is the next angular momentum value that the system can take. In our calculations 
we choose to use the Harris formula to the third order. 

At this point we have to specify what form the coupling will take. The exact form of 
the inter-band interaction is not known. We have chosen to perform the calculations using 
the standard band coupling formalism 0. Since the ground state has K = 0, the possible 
cases are limited by the values of K chosen for the excited band. In the following we allow 
K = 0A and 2. 



For K = ^ K = 1 
For K = 0^ K = 
Constant interaction 
For K = 0^ K = 2 



C{I) = C^I{I+1), (9) 

CiI) = CI{I+l), (10) 

= C. (11) 

C(/) = C(/ + l)(/ + 2)(/-l). (12) 



The actual angular momentum dependence of the coupling term is not decisive for the 
following reason: The mixture of two uncoupled band states with the consequent distortion 
of the shapes of the energies as functions of the angular momentum depends on the strength 
and the form of the coupling term as well as on the energy difference of the unperturbed 
states. The bands will be close to each other (in the Energy vrs. J plot) only for a single 
value of angular momentum, and therefore the mixture will be insignificant for all other 
values of angular momentum. As soon as the strength of the coupling term is fitted at 
this angular momentum, the actual functional dependence of the coupling strength becomes 
irrelevant. For the actual calculations we tried different forms of the coupling and we state 
in Table | which ones yield the best fit for the nucleus under study. 

At this point we have a complete phenomenological description for the backbending 
phenomena. If we use to the expansion Eq. (P) up to the third order, we then have eight 
free parameters to fit, namely the six coefficients for the two uncoupled bands, ^'j', A\, Al, 
A2, Al and A2, the band-head energy Eq, and finally the coupling C. This is not a difficult 
task, and we will describe the details and show the results in the following section. 

To calculate the BE{2)^s we have to express the matrix elements of the quadrupole 
operator as a function of the matrix elements in the uncoupled system. To be more specific 
with the notation, states of the unperturbated Hamiltonian that belong to the ground-state 
band will be denoted by |/M0) and states that belong the excited band \IMK). For the 
coupled system states have good angular momentum /, its projection M, but K is not a good 
quantum number any more. We denote states of the complete system as \IM1) and \IM2). 
After diagonalizing the Hamiltonian, we can express the eigenstates of the full system as a 
linear combination of the unperturbed states 

|/M1) = Ao(/)|/MO) + Xi{I)\IMK), 

\IM2) = /io(/)|/MO) + fii{I)\IMK), (13) 
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where Ao(/), Ai(/), fioil) and will be provided from the diagonalization of the Hamil- 

tonian. Consequently we can express the matrix elements of the quadrupole operator, in 
the coupled function of the matrix elements of the unperturbated case, 



IMl\Q\I'M'l) = Ao(/)Ao(/') [IMO\Q\I'M'0) + Ai(/)Ai(/') {I M K\Q\I' M' K'^ 

+ Ao(/)Ai(/') (lMQ\Q\rM'K) + Ai(/)Ao(/') {l M K\Q\I' M' 

IM1\Q^\I'M'2) = Ao(/)/io(/') {lM0\Q\I'M'0) + Ai(/)/ii(/') (l M K\Q\I' M' k) 

+ Ao(/);Ui(/') {iMQ\Q\I'M'k) + Ai(/)^o(/') {lMK\Q\I'M'{]) , 

IM2\Q^\I'M'2) = /io(/)/io(/') (/M0|Q|J'M'0) + (/Mir|Q|/'M'ir) 

+ /io(/)/xi(/') (/MO|Q|/'M'ir) + /ii(/)/io(/') (/Mir|Q|/'M'o; 

(14) 

The reduced matrix elements of the quadrupole operator Q for the uncoupled case are taken 
from the phenomenological Bohr-Mottelson model. 



JO II Q II I'k) = aK> g'^^'^'^^-^^'^°^2y27TTy27^ 1^'^ ) , (15) 




aK' 




K' = 



where gbandi^band2 ^^^^ ^^^q intrinsic quadrupole moments of any inter or intra-band tran- 
sitions. Consequently we have three more parameters to fit, namely the three intrinsic 
quadrupole moments. From measured BE{2)^s we fix these three parameters and then 
calculate the rest of the BE{2)^s that are needed in the calculation but not available exper- 
imentally. 

For the first application, ^^^W, we need the excitation energies and the quadrupole matrix 
elements of the neighboring even-even nuclei ^'^^'isow. Since we assume particle-number non- 
conservation the difference between the two neighboring nuclei should be small, and since 
there are not enough experimental data for ^''^W, we did the calculations for ^'^'^W only. In 
the special case where the excited band has odd values of angular momentum the fitting can 
be simplified further. These odd angular momentum states do not interact with the states 
of the ground state band, and since they belong to the same band, the parameters A of 
this band can be calculated independently. First, we fit the parameters ^exc^ the states 
with odd value angular momentum, and then we fit the other 5 parameters. For example 
the excited band of ^^°W (see Fig. (|1|)) has states with odd value /. Therefore, we fitted the 
upper band first and when the ^eic^ were found, we fitted the rest. The results are shown 
in Table (p along with Fig. (|I|). 

For the BE{2)^s we utilize the few experimental data available. As stated above, there 
are three parameters to fit, go, Qk and qko- With these parameters fixed, we calculated all 
transitions that will be needed for the calculations of the odd nuclei. The results of the fit 
are given in Fig. (Q) and Table (|). 

As a second application we considered the odd proton nucleus ^^^Lu. The neighboring 
even nuclei are ^^^Yb and ^^^Hf. There are enough experimental values for both these 
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nuclei, and therefore we have fitted both nuclei and taken the average values of the excitation 
energies and quadrupole matrix elements for application to the odd nucleus. The calculation 
was done as described above for ^^°W with one modification however. We realized that the 
fitting of the BE{2) values is more sensitive to the strength of the coupling C{I) than are 
the excitation energies. Therefore, we first do a preliminary fit of the excitations energies. 
We then adjust the wavefunctions found from these calculations to fit the BE{2) values by 
allowing C to vary once more, in addition to the necessary choice of intrinsic quadrupole 
moments. The results are shown in Fig. (0) and in Table (1). Then the excitation energies 
are re-calculated with an already fixed coupling strength C, and we find new values for 
the ^'s. The results of final fitting of the excitation energies are shown in Fig. (^) and in 
Table (D. 

III. CALCULATIONS AND RESULTS FOR ^^^W 

The equations of motion (EOM), as well as the normalization conditions, for the case of 
backbending are not formally different from those presented in the two previous applications 
lU and therefore will not be repeated. In contrast to these applications, the quadrupole field 
F will not take one of the simplified forms of the Bohr-Mottelson theory, such as Eq. (P^. 
Instead, it will take values that are determined from the calculation described in the previous 
section, culminating in Eq. (Q). The form of F is given by, 

Tain,crn';j = -\ ^ J | (/n||g||/V) , (16) 

where n = 1, 2 , ra' = 1,2, and the quadrupole matrix elements, {Jn\Q\l'ri!) ^ were calculated 
in the previous section. Notice that F depends on the total angular momentum J in contrast 
to the two previous cases. In the previous calculations, because we assumed axial symmetry, 
F was independent of J. In this application the symmetry is lost, with implications that will 
be discussed later. 

Similarly to the previous two examples, we used a large single-particle space (including 
all states from 5 major shells). The energies and matrix elements of these single-particle 
levels were calculated using the Woods-Saxon potential. The core excitations oji were taken 
from calculations for the core described in the previous section. At this point all input 
parameters are fixed and the only thing remaining to be done is to solve the eigenvalue 
problem given by the EOM. In the same way as before, the strength of the quadrupole 
field is treated as a free parameter and the values of the single-particle energies found from 
Woods-Saxon calculations are allowed to vary by ±5%. 

Two remaining technical problems should be discussed before we present the results 
of the calculations. The first difficulty in solving the EOM is that the set of solutions is 
over-complete by a factor of two, as has previously been discussed. This is a consequence 
of the fact that the basis states form an over-complete (and, therefore, non-orthogonal 
set). Consequently half of the states found by solving the EOM are not physical and have 
to be removed. A technique has been developed to do so and was presented in [||]. In 
0, an extension of the technique to the case of avoided crossings, was discussed. In this 
paper, we summarize only the essential points. The Hamiltonian is first decomposed into 
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anti-symmetric and symmetric parts with respect to particle-hole conjugation. If only the 
anti-symmetric part is diagonalized, then for every positive energy eigenvalue there is a 
negative partner. From the BCS theory we know that the positive eigenvalues are the 
physical solutions and the negative solutions the non-physical ones. Then the symmetric 
part was turned on "slowly" and at every step the physical solutions were identified using 
the projection operator for the physical space built from the wavefunctions of the previous 
step. The projection operator has an eigenvalue 1 for all real states and for the non- 
physical states. Since this operator is constructed from wavefunctions of the previous step 
we expect that it will give an expectation value close to 1 for physical and close to for the 
non-physical states of the current step. This procedure continues until the symmetric part 
is turned on fully. In we explained the case where physical and non-physical states of the 
same angular momentum come very close to each other, and consequently the corresponding 
wavefunctions change rapidly. In this instance the projection operator method is not valid 
unless the steps taken are extremely small, making the actual calculation numerically very 
slow. In this case we utilized the fact that states of the same angular momentum can not 
actually cross. As stated previously, the two bands involving the near crossing interchange 
the characters of their wavefunctions as they pass near each other. Since the projection 
operator method is basically a comparison of wavefunctions, the form of this operator just 
prior to the crossing would classify a physical state as unphysical and vice versa when 
used in coarse steps involving the near-crossing. Therefore when this occurs, we adjust 
the projection operator appropriately in order to make the correct identification in future 
steps. The problem in now reduced to determining if there is a crossing. This was done by 
comparing the energy differences between two states at two consecutive steps. This method 
worked extremely well both in the example of the previous work and in the present 
case. 

The second technical problem is the classification of states into different bands. In the 
first two applications where K was a good quantum number, the identification of bands was 
done based on the K value of the band. The technique was explained in [Q] and again we 
repeat only essential aspects. In the case that the core is approximated as a rigid rotor, 
and the formulas of the geometrical model are valid, it turns out that the quadrupole field 
is independent of the total angular momentum J. This is because of the assumed axial 
symmetry. As a result when only the antisymmetric part of the Hamiltonian is considered 
the Hamiltonian has no J dependence. It follows that for a given K, i.e. a given band, 
all eigenstates with different J are degenerate. We first solve our equation for a minimum 
possible angular momentum, J = \- This identifies the K = ]^ bands. At J = |, we find a 
set of energies equal to those previously found and possibly additional solutions identified 
as the band heads for K = ^. Continuing in this way we assemble all solutions into a series 
of fiat bands with identified K values. 

In the case that backbending occurs, even when the core excitations are neglected (sym- 
metric part of H set to zero) the remaining Hamiltonian is J dependent. Therefore, we can 
not use the method described above. We can however, take the additional steps of turning 
off the coupling between the two core bands to return to the axially symmetric case. Then 
the previous procedure can be used, but there is the additional complication of having to 
turn on both the symmetric part and the coupling between the two core bands. We can do 
that in any order, i.e. we can first turn on the coupling between the core bands and then the 
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core excitation (symmetric part of the Hamiltonian) or the other way around. We choose 
to do the former. First we turn on the band couphng slowly, and as a result we break the 
axial symmetry. At every step we compare the wavefunctions at two consecutive steps by 
calculating the overlap integral. It turns out that even though the J degeneracy is lifted the 
wavefunctions of the axial and non-axial cases (not including the core excitations) are very 
similar. Only for a few eigenstates was the mixture between states belonging to different K 
bands so big we had to turn on the coupling slowly. In most cases large steps were sufficient. 

The first application made was to the nucleus ^''^W. Recent observations [lTO| , pl| have been 
interpreted as showing that for this odd-neutron nucleus, an alignment to an axis intermedi- 
ate between the deformation axis and the collective rotational axis (Fermi alignment [Q), 
gives rise to a backbending at J ~ (Frauendorf has called the associated rotational 
structure a t-band, since the cranking- mo del description requires a tilting of the cranking 
axis away from the principal axes of the intrinsic prolate spheroidal shape.) Nevertheless we 
have chosen to try the standard two-band model. 

Energy levels of ^^^W calculated from the present work are presented in Fig. (^. As 
can be seen from the figure, the yrast band {K = 7/2) is reproduced with high accuracy. 
The same can be said of the first excited band {K = 7/2 or tilted-band). Then follow two 
K = 1/2 bands which agree very well with the theory. The most striking feature is the exact 
reproduction of the staggering nature of the first K = 1/2 band. At about 1 MeV above 
the ground state, theory predicts a = 5/2 band which experimentally is not observed. 
We believe that this a weakness of the fitting routine for the following reasons. The relative 
band-head energies were fitted, using a standard fitting routine, to the observed values by 
varying the strength of the quadrupole field, k. The fitting was done by minimizing the x^- 
In the absence of an experimental value and the presence of a a theoretical prediction the 
contribution to the is 0. In other words for the fitting routine (at least for the one we 
used), if the theory predicts a band at low energy that is not observed experimentally it 
does not increase the x^- We are currently working to build a better fitting routine which 
will be able to avoid such a problem. The result shown in Fig. ^ already represents some 
improvement over the fitting procedure used initially. 



IV. RESULTS FOR ^^^LU 

The second backbending application we tried was to the proton spectra of ^^^Lu. This 
nucleus was studied previously by methods related to those of this paper by Chen et al []T3| . 



Since it is generally accepted that the s-band in an even nucleus is formed by the decoupling 
of a pair of quasiparticles from the ground-state band, these authors, among others, decided 
that it was necessary to use a semi-microscopic description of the backbending in the neigh- 
boring even nuclei. However, their fit to the spectra above the backbend leaves something 
to desired. For this reason, we have chosen to separate the problem of fitting the backbend 
in the even nuclei from that of fitting the corresponding spectra in the odd nuclei by using a 
phenomenological two-band mixing model for the even nuclei, which allows, simultaneously 
an accurate fit to the the spectrum of the crossing bands and to all observed B(E2) values 
(see previous section). In Fig. (|^), we compare our fit for both negative and positive signa- 
ture (a = — |, a = |) crossed bands in ^^^Lu with those of Chen. The closeness of our fit 
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compared to that of the previous authors appears simply to reflect corresponding deviations 
in the spectra of the even neighbors. 

We next study magnetic transitions of the yrast band and compare the result to experi- 
ment and other theories. In we have presented the equations for the Ml matrix element 
and calculations for the B{M1) transitions of the ground state band of ^^^'^^''Gd. To calcu- 
late B{M1) rates, we must first evaluate the matrix element of the Ml operator between 
states of the even neighboring cores (/||M||/) and the matrix elements of the same operator 
between single particle states {rriac)- In the case of (/||M||/), and because of the limita- 
tion on the availability of the experimental data, we have chosen to use a phenomenological 
model for the value of g, that is derived in [0], namely 



9 = 9r+ {9j - 9R)i/I, (17) 

where gn refers to the ground band, i is the aligned angular momentum, and gj is the single 
particle g factor for the j shell in which alignment occurs. For the 1^13/2 level, gj takes 
the value of —0.20 and gn = 0.31. In Fig. ^ we show the g factors for the yrast band of 
^^^Hf obtained from the cranking model, the two-quasiparticle plus rotor model, and from 
the simple formula shown above. For the single particle matrix elements we used a simple 
formula from Eq. (44) of ||]. 

Results from different calculations of 5 (Ml) rates of the yrast band of ^^^Lu are summa- 
rized in Fig. ^ Figure |] shows experimental values from [0, the result from the cranking 
calculation of [jl3l, a core-particle coupling model also presented in [|l3l and finally our cal- 



culations. Our work reproduces the experimental values quite well and clearly better than 
the other works. 



V. CONCLUSIONS 

In this paper we have presented yet another distinct application of the KKDF model. 
Calculations for ^^^Lu and ^^^W were performed starting from a phenomenological descrip- 
tion of the neighboring even nuclei. The results are sufficiently satisfactory to encourage us 
to consider still more sophisticated applications of the model utilized. A fundamental ap- 
proximation that all applications so far shared was that the incorporation of particle-number 
conservation into the equations has little effect on the results, because the properties of the 
neighboring nuclei are so smooth and slowly-varying with particle number. It therefore be- 
hooves us to study a case where the incorporation of exact number conservation is essential 
because we are in a transitional region. Such an has been worked out W7 . 
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TABLES 



nucleus 


-4^g 






A^ 


A^ 

•^exc 


•^exc 


Band-Head 






r 1 




\-keV\ 


r 1 


r ?i6 1 


[keV] 


180 1^ 


17.16 


-2.19 10-2 


1.95 10-5 


6.45 


4.00 10-3 


1.7 10-6 


1663.50 


164^5 


18.79 


-2.56 10-2 


2.10 10-5 


6.03 


2.21 10-3 


2.10 10-5 


1603.34 




20.55 


-2.62 10-2 


1.90 10-^ 


6.10 


-2.62 10^ 


1.92 10-' 


1675.34 



nucleus 


Qo 


QK 


QOK 


Coupling Type 


Coupling 








m 


[keV] 


[keV] 


180^ 


0.711 


0.031 


0.66 


C(/)=C(/ + l)(/ + 2)(/-l) 


4.21 10-5 




0.540 


0.021 


0.41 


C{I) = CI{I+1) 


1.30 10-1 




0.560 


0.019 


0.39 


C{I)=CI{I + 1) 


1.37 10-1 



TABLE I. Values obtained for the parameters for the phenomenological model of the even 
cores. The subscribe "exc" on Aexc identifies the values of the excited bands. 
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FIG. 1. Fit of the ground-state band and the first excited band of ^^°W to a phenomenological 
two-band model. Solid lines represent the theoretical values with coupling equal to 4.3 10~^keV. 
The experimental values are represented by circles. 
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FIG. 2. Experimental and theoretical BE(2) values for the yrast band and first excited band. 
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FIG. 3. Energy levels of ^^^Yb (left) and ^^^Hf (right). Both experiment and theory are shown. 
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FIG. 4. BE{2) values from experiment (circles) and theoretical fit (solid line) of ^^^Yb. 
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FIG. 5. Comparison of theoretical results to observed energy levels of ^^^W. 
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FIG. 6. Comparison of the observed yrast band of Lu with theory. On the left are the 
negative signature states (a = — i) and on the right are the positive signature states (a = i). 
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FIG. 8. B{ml) transitions rates in /i^. On the top are the positive signature states {a = —\) 
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